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THE QUASAR SDSS J1536+0441: AN UNUSUAL DOUBLE-PEAKED EMITTER 
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ABSTRACT 

The quasar SDSS J153636. 22+044127.0, exhibiting pecuhar broad emission-hne profiles with mul- 
tiple components, was proposed as a candidate sub-parsec binary supermassive black hole system. 
More recently, imaging revealed two spatially distinct sources, leading some to suggest the system to 
be a quasar pair separated by ~5 kpc. We present Palomar and Keck optical spectra of this system 
from which we identify a third velocity component to the emission lines. We argue that the system is 
more likely an unusual member of the class of active galactic nuclei (AGNs) known as "double-peaked 
emitters" than a sub-parsec black hole binary or quasar pair. We find no significant velocity evolution 
of the two main peaks over the course of 0.95 yr, with a 3ct upper limit on any secular change of 
70 km s~^ yr~^. We also find that the three velocity components of the emission lines are spatially 
coincident to within 0'.'015 along the slit, apparently ruling out the double-quasar hypothesis. 
Subject headings: quasars: individual (SDSS J153636. 22+044127.0) — accretion, accretion disks 



1. INTRODUCTION 

When galaxies merge, it is expected that their cen- 
tral supermassive black holes will form a binary which 
will also eventuall y merge via the em ission of gravita- 
tional radiation (Bcgcl man et al.lll98 0l. The process by 
which the orbit of the two black holes shrinks to the 
point at which gravitational radiation is effective in pro- 
ducing a merger within a Hubble time is not currently 
understood (dubbed the "fi nal parsec problem"; e.g., 
iMilosavljevic fc Merrittl[2003l ). and the observational sig- 
natures of such systems have proved to be elusive. 

Boroson & Lauer (2009; hereafter BL09) recently iden- 
tified a candidate sub-parsec binary supermassive black 
hole system which could potentially merge in less than 
a Hubble time. They searched the S loan Digital Sky 
burvey (SDSS; Abazaiian et al.l [20091 ) quasar spectra 
and identified the object SDSS J153636.33+044127.0 
(hereafter SDSS J1536+0441) as one of two low-redshift 
{z < 0.7) quasars exhibiting components having multi- 
ple redshifts. The key features of the system are two 
components of the broad emission lines separated by 
~ 3500 km s^^. They interpreted these to be the result 
of broad-line regions around a pair of black holes with 
masses of lO'^'^ and 10®'^ M0, with the relative velocity 
being a result of the orbital motion in a binary having a 
sub-parsec semimajor axis with a period of ^^100 years. 

The night after the BL09 publication came to our at- 
tention, we obtained Palomar spectra and identified a 
third velocity compon ent of the broad emission lines 
(jChornock et al.l l2009f ) . This led us to conclude that 
SDSS J1536+0441 was instead an unusual member of 
the class of AGNs known as "double-peaked emitters" 
(e.g., Halpern & Filippenko 1988; Eracleous & Halpern 
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1994; Strateva et al. 2 003), a suggest ion that was later 
made independently by Gaskell (2009[). 

Radio observations by Wrobcl & Laoil (j2009f ) led those 
authors to suggest a third possibility: that of a quasar 
pair. They identified two radio sources separated 
by 0'.'97 in Very Large Array ( VLA) observa t ions o f 
SDSS J1536+0441. Subseauentlv. iDecarh et all (|2009aD 
found that the two radio sources were coincident with a 
pair of K-hand sources that they suggested were mem- 
bers of a '^5 kpc quasar pair. As we completed this work, 
a second preprint by Lauer & Boroson (2009; hereafter 
LB09) appeared. They presented Hubble Space Telescope 
(HST) images confirming the double nature of the opti- 
cal source, but they also argued that the optical data 
were inconsistent with a double-quasar hypothesis. 

In §2 of this Letter, we present the Palomar spectra 
that led to our initial identification of this source as a 
double-peaked emitter, as well as additional high-quality 
Keck spectra. We set upper limits on the long-term evo- 
lution of velocity shifts of the broad peaks and on any 
spatial offset of the emission regions in §3. We then ar- 
gue (§4) that the double-peaked emitter nature of this 
object provides a better explanation than the alterna- 
tives. LB09 reached similar conclusions, so we compare 
with their results where appropriate and highlight our 
differences. 

2. OBSERVATIONS 

Our initial spectroscopy was performed using the Dou- 
ble Spectrograph o n the Palomar 5-m Hale telescope 
(|Oke fc Gunnlll98l) : two 600 s exposures were obtained 
on 2009 March 7.52 UT. The blue-side (3150-5700 A) 
observations we present here had a spectral resolution of 
'^B A from the combination of a 1'.'5 slit and the 300/3990 
grating. In addition, we obtained a pair of 600 s observa- 
tions with the Echellette Spectrograph and Imager (ESI; 
Sheinis et al. 2002) on the Keck II 10-m telescope on 
2009 March 22.64 UT. The 0'.'75 slit was used, giving a 
resolution of ~45 km s~ over the range 4000-10,200 A. 

The two-dimensional (2D) images were processed and 
the spectra were extracted using standard routines in 
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IRAF^ . Our own IDL routines were used to flux calibrate 
the individual ESI echellette orders and remove telluric 
absorption features, as well as to combine the individ- 
ual orders together over their overlap regions to make 
the final spectrum. The wavelength scale was estab- 
lished using comparison-lamp spectra, with small shifts 
(< 0.05 A) of the zeropoints of the wavelength scales 
found by cross correlating the night-sky lines in the ob- 
ject observations against reference-sky spectra. The flux 
level of the first ESI observation was scaled up by 15% 
to match the second before combining them, under the 
assumption that the first was affected by clouds or poor 
guiding. 

3. RESULTS 

Both of our spectra reveal that SDSS J1536-H0441 has 
unusual broad-line profiles (see Fig. [T]), as initially de- 
scribed by BL09. The Balmer lines have a multi-peaked 
structure, with one peak having the same redshift as the 
narrow nebular emission lines {z = 0.3889) and a second 
one blueshifted by ^ 3500 km s~^ relative to the first. 
We follow BL09 by referring to these as the "r-system" 
and "b-system," respe ctively. However, as we first noted 
(jChornock et al.ll2009[ ). the Ha line shows a clear third 
component, visible as a shoulder or hump present at a 
velocity of ~4000 km s~^ relative to the r-system, al- 
most symmetrically opposite in velocity space from the 
b-system. The third component is also present in the 
H/3 profile, but it is partially hidden by the superposed 
[O III] A4959 emission line. Although the broad Mg II 
A2800 line is contaminated by superposed narrow Mg I 
and Mg II absorption lines as well as broad Fe II emis- 
sion, it too shows a boxy profile with the "corners" of the 
box at the same velocity as the b-system and the third 
component. 

No features due to starlight in the host galaxy are seen 
in the spectra, but absorption lines from several com- 
mon interstellar species (Na I, Ca II, Mg I, and Mg II) 
are present at a redshift intermediate between the b and 
r-systems (z = 0.3878), as identified by BL09. Our ESI 
spectrum, which includes only the Na I and Ca II lines, 
resolves those absorptions (full widths at half maximum 
w 100 km s""'^ and 170 km s""'^, respectively). Weaker 
(rest- frame equivalent width w 0.16 A) absorption from 
Na I D2 is also present aX z — 0.3890, close to the redshift 
of the r-system, and is likely due to the host galaxy (Dl 
at this redshift is blended with D2 of the stronger absorp- 
tion system). The stronger absorption system is there- 
fore blueshifted by 220 km s~^ relative to the narrow 
emission lines, possibly analogous to the blueshifted low- 
ionization ultraviolet (UV) absorber see n in the proto- 
typical double-peaked emit ter, Arp 102B (jHalpern et al.l 
ll996HEracleous et al.l[200l . 

It is of interest to examine the line profiles of 
SDSS J1536-I-0441 in more detail, both to investigate po- 
tential line-profile variations and velocity shifts of emis- 
sion components. In Figure [21 we have plotted a com- 
parison of the Balmer-line profiles in the SDSS and ESI 
spectra. The original SDSS spectrum did not cover the 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1.— Emission-Une profiles of SDSS J1536-I-0441. Shown are 
Ha (black) and H/3 (blue) from the Keck/ESI spectrum and Mg II 
A2800 (red) from the Palomar spectrum. The moderate-resolution 
ESI data have been rebinned for display purposes. A simple spline 
fit to the continuum beneath each emission line has been removed. 
The zeropoint of the velocity scale was set by the narrow [O III] 
line associated with the r-system. Superposed narrow emission and 
absorption lines from other species are shown as dotted lines. 

full Ha profile, explaining why BL09 could not have iden- 
tified the third component. The spectra are very simi- 
lar, with the biggest difference being the ~10% higher 
peaks associated with the b-system in the newer spec- 
trum. This increase in the ratio of the peak fluxes of the 
b-system to the r-system is seen consistently in Ha, H/3, 
and H7. 

The lower- left panel of Figure [2] also shows a compari- 
son of the blue and red halves of the Ha profile in our ESI 
spectrum. There is perhaps some excess flux at high ve- 
locities on the blue side relative to the red half of the line, 
but the overall shape of the wings {\v\ > 5000 km s^^) is 
similar on both sides. The big difference between the two 
halves is the emission peak at —3500 km s~^ associated 
with the b-system, which is at a very similar velocity (but 
with the opposite sign) as the shoulder on the red side. A 
potential explanation is that the b-system and the third 
component represent velocity flows in the broad-line re- 
gion that are projected onto our line of sight in opposite 
directions (possibly due to rotation or outflows). The 
binary black hole interpretation lacks a natural explana- 
tion for both the existence of the third component and 
its presence at an opposite velocity from the b-system. 

One explicit prediction of the binary black hole hy- 
pothesis made by BL09 is that the system might show 
noticeable velocity shifts due to the relative orbital mo- 
tion of the two components of order 100 km s^^ yr^^. 
Our ESI spectrum was taken 0.95 yr after the SDSS spec- 
trum, so such a velocity shift should be easily measur- 
able. We experimented with several methods to find a 
robust method of determining line peaks without being 
sensitive to assumptions about the line profile. We set- 
tled on a procedure where we iteratively fit a quadratic 
function to each line peak and then fit a quadratic to 
the region of ±500 km s^^ around the preliminary peak 
location. The two spectra were fit in an identical man- 
ner. The differences between the peak wavelengths in 
the SDSS and ESI spectra were then converted to veloc- 
ity shifts and are listed in Table [H Gaussian fits to the 
peaks gave similar results for the shifts. 

The two spectra have very consistent wavelength 
scales, as demonstrated by the 2.8 km s^^ shift between 
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Fig. 2. — Comparison of Balmer-line profiles. Our Keck/ESI 
spectrum of SDSS J1536-I-0441 is shown in black and the SDSS 
spectrum analyzed by BL09 is in red. The lower-left panel com- 
pares the blue portion of the Ha profile in the ESI spectrum (black) 
to the red half, reflected across zero velocity (red-dashed). 



TABLE 1 

Offsets between SDSS and ESI Spectra 



Line 



Velocity Shift (km s^^) 



[O III] A5007 r-system 


-2.8 


± 


2.5 


H/3 r-system 


13 


± 


18 


Ha r-system 


4.8 


± 


27 


H/3 b-system 


4.8 


± 


22 


Ha b-system 


-19 


± 


15 



the [O III] lines in the two spectra. Our ESI spectrum 
has a high signal-to- noise ratio (S/N) and velocity resolu- 
tion, so the errors are dominate<i by the S/N of the SDSS 
spectrum, as demonstrated by the substantially smaller 
errors in fitting the bright [O III] line. The weighted 
average values for the velocity change of the b and r- 
systems (as measured from the Ho; and H/? lines) are 
—18 ± 19 km s"^ yr~^ and 16 ± 23 km s~^ yr~^, respec- 
tively, in the rest frame of the object. We can therefore 
set a 3ct limit on any velocity shift of the underlying 
e mission regions o f 5,70 km yr~^. 

iWrobel fc Laml (|2009l ) observed SDSS J1536-H0441 
with the VLA and discovered a pair of radio sources 
with 0'.'97 separation (named A and B). Decarli et al. 
(2009a, b) and LB09 found if-band and optical sources 
coincident with both the A and B radio components. 
One possible interpretation of these results is that the 
A and B radio sources can be identified with the b and 
r-systems seen in the optical spectra, making this source 
a double quasar with a ^5 kpc separation. With a fiber 
size of 3" in diameter, both systems would have naturally 
fallen within the single spectrum obtained by SDSS. At 
the time of our E SI observations, we w ere unaware of 
the radio results of lWrobel fc Laoij (|2009[ ) and so did not 
orient the spectrographic slit along the vector separat- 
ing the A-fB radio sources (PA=88°). Instead, our slit 
was aligned with the parallactic angle (Filippenko 1982; 
PA=48°).We also do not know for certain the exact posi- 
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Fig. 3.— Spatial profiles of SDSS J1536-I-0441 at several difll'erent 
wavelengths. The profiles for the order containing H/3 (left panel) 
represent sums of the counts over wavelength intervals correspond- 
ing to the peaks of the b and r-systems, the [O III] A5007 line, as 
well as continuum regions on either side of the emission complex. 
The spatial profiles for the order containing Ha (right panel) rep- 
resent the peaks of the three components labeled in Figure [l] as 
well as a continuum region blueward of the emission-line complex. 

tion of the slit on the sky relative to the two components. 

Whether the slit was centered on the bright optical 
source corresponding to the A source or somewhere be- 
tween the A and B sources, one would expect an easily 
measurable spatial offset between the b and r-systems 
of 0'.'7 projected along the spatial direction of the slit 
if they were associated with the A and B sources. The 
echelle orders of ESI are highly curved relative to the 
CCD columns, so we used the trace of our bright stan- 
dard star to rectify each order. We summed the counts in 
the 2D CCD frame in several wavelength intervals corre- 
sponding to the three velocity components of the broad 
emission lines, as well as the narrow [O III] A5007 line 
and the surrounding continuum. The results for the ESI 
orders containing Ho; and II/3 are plotted in Figure [3l 

As can be seen from Figure |3l the spatial profiles are 
almost identical at all wavelengths along the slit. We fit 
Gaussians to the spatial profiles and found the centroids 
to be essentially coincident. The maximum and mini- 
mum slit positions across the Ha and H/3 line profiles 
are equal to better than 0.1 ESI spatial pixel (< 0'.'015), 
which is highly inconsistent with a 0'.'7 projected separa- 
tion. Therefore, we reject the double-quasar hypothesis. 
LB09 obtained a similar constraint on the basis of obser- 
vations taken with the slit aligned along a more optimal 
position angle and also presented photometric evidence 
that the B source is associated with an early-type galaxy. 

4. DISCUSSION 

AGNs with double-peaked broad emission lines have 
been known for many years as a relatively rare sub- 
class of AGNs. Systematic surveys of radio galaxies 
(lEracleous fc H alpcrn 1994) and the SDSS quasar sam- 
ple (jStrateva et al...2003 ) have found numerous members 
of the class. The most plausible current explanation for 
the displaced red and blue peaks for most objects in this 
class inv olves rotational m otion in a relativistic accre- 
tion disk (|Chen et al.lll989l ). possibly with some emission 
component originat ing in a wind being driv en from the 
surface of the disk (|]V[urray fc Chian3ll997f ). 

Some double-peaked emitters observed in the past have 
been proposed to be candidate binary black hole sys- 
tems (|Gaskelll ll983. 1996), much like SDSS J1536+0441. 
However, that interpretation has been ruled out in well- 
studied systems such as Arp 102B, 3C 390.3, and 3C 332 
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by the lack of appropriate long-term velocity evolution of 
the p eaks (jHalpern &: Filippenkd 119881 : lEracleous et al.1 
|1997() . BL09 proposed that the observable velocity evo- 
lution would be a test of their black hole binary hypoth- 
esis. While the observed upper limits (§3; LB09) to the 
velocity evolution constrain the allowed orbital parame- 
ters for a potential binary (explored in detail by LB09), 
we caution that small velocity shifts of the line peaks may 
even be expected in th e double-peaked e mitte r scenario. 
The objects studied bv lEracleous et al\ (jl997) did show 
significant velocity evolution, but in a manner that was 
inconsistent with orbital motion. T he origin of these va ri- 
ations is still not well understood (jGezari "eral][200l . 

The full widths at half maximum and at quarter max- 
imum for Ha in SDSS J1536+0441 (relative to the max- 
imum flux in the central peak) are 8700 and 11,500 
km s^^, respectively. These values are almost identical 
to the mean values o f 8852 and 11,702 km s~^ found by 
IStrateva et al.l (|2003l ) for their optically selected double- 
peaked emitter sample from SDSS, and are significantly 
larger than those of their reference sample of normal 
quasars. 

The [O III] A5007 flux in the ESI spectrum is within 
2% of that in the SDSS spectrum, showing excel- 
lent agreement of the absolute flux scales despite the 
very different aperture sizes (the A5007 emission re- 
gion appears to be spatially unresolved; Fig. [3|). The 
total continuum-subtracted Ha flux over the velocity 
range \v\ < 10,000 km s~^ is 1.6 x lO^-^'^ ergs cm~^ 
s~^, or Liia — 8.2 X lO''^ ergs s~^ (assuming a flat 
ACDM cosmology with Hq = 72 km s"^ Mpc"^; 
Dunkley et al. 2009). This value is higher than 
any of the double-pe aked emitters from the sample of 
IStrateva et al.l ()2003f ). but is exceeded by a few of the 
radio-selected objects of lEracleous fc HalpernI (j2003D . in- 
dicative of the extreme nature of SDSS J1536-I-0441. In 
addition, the low-ionization narrow emission lines have 
lower rest-frame equivalent widths (~0.7 A for [O I] 
A6300 ) than most double-peaked emitters in the lite r- 
ature (lEracleous fc Halpernll200llStrateva et al.ll2003f l. 

LB09 compared SDSS J1536-h0441 with several ob- 
jects from the literature, identifying three differences be- 
tween SDSS J1536-f0441 and the majority of double- 
peaked emitters. Two of these differences are the rel- 
ative strengths and shapes of the central and blue peaks 
seen in SDSS J1536-I-0441 compared to the other objects. 
Doubl e-peaked emitters are rather heterogeneous as a 
class (jEracleous fc HalpernI fl99l IStrateva et al.ll2003f ). 
and we suggest that SDSS J1536-I-0441 is an extreme ex- 
ample. The third difference noted by LB09 was the lack 
of change in the intensity of the blue peak over a one-year 
baseline seen in their Kitt Peak spectrum. However, in 
our data (Fig. [2]) the blue peak has increased in flux by 
~10% in both Ha and H/3. 

5. CONCLUSIONS 

We demonstrated in §3 that the broad emission com- 
ponents of SDSS J1536-I-0441 were spatially coincident, 
ruling out the double -quas ar hypothesis suggest ed by 
iWrobel fc Laa^ (l2009h and iDecarh et al] (|2009aD . We 
believe the preponderance of evidence is in favor of 
interpreting SDSS J1536-I-0441 as an unusual double- 
peaked emitter, not the binary black hole originally sug- 
gested by BL09. LB09 have reached similar conclusions 



about the nature of this system, but do not yet rule out 
the binary hypothesis. As a phenomenological matter, 
SDSS J1536-I-0441 clearly shows multi-peaked emission- 
line profiles like other members of the class of double- 
peaked emitters. We cannot conclusively disprove the 
binary black hole hypothesis, but feel that there is no 
compelling reason to adopt a binary interpretation when 
SDSS J1536-I-0441 shows characteristics of a previously 
known class of objects probably unrelated to binaries. 

SDSS J1536+0441 certainly exhibits a few extreme 
properties for a double-peaked emitter. The optical lumi- 
nosity is very high as is the ratio of the blue to red peaks. 
iTang fc Grindlayl ()2009l ) modeled the Balmer line profiles 
in the Keck data presented here and found that a stan- 
dard disk model could reproduce the third component 
and a shelf of emission under the peak associated with 
the b-system. However, their disk models alone could 
not reproduce the sharp blue peak, so another compo- 
nent was required, which they proposed was due to a 
broad-line region around a secondary black hole in the 
accretion disk. We note that a significant fraction of the 
other known double-peaked emitters are also not well fit 
in detail by simple disk models, either, possibly indicat- 
ing the presence of an add itional emission component 
or non -axisymmetric disks (lEracleous fc HalpernI 120031 : 
IStratev a et al. 2003). 

Instead of waiting for several years to further refine 
the constraints on the lack of orbital motion in the sys- 
tem as suggested by LB09, we note that there is poten- 
tially a simple test that can be done more quickly. Now 
that HST has been refurbished, it will be possible to ob- 
tain UV spectroscopy. UV observations of other double- 
peaked emitters have shown that the emission lines from 
high-ionization material (e.g., Ly a and C IV) frequently 
lack double-peaked velocity profiles ("Halp ern et al.in~996l : 
lEracleous et al. 2004). This striking difference between 
the profiles of different lines has been attributed to den- 
sity and radiative-transfer effects in the multiple line- 
emitting regions (e.g., the accretion disk and associated 
wind). The absence of multiple velocity peaks in the 
UV emission lines would be strong evidence in favor of 
SDSS J1536-I-0441 being a member of the class of double- 
peaked emitters; the presence of multiple peaks, how- 
ever, would still render the interpretation ambiguous. 
UV spectra may also clarify the nature of the strong 
blueshifted absorption system. 

The high spatial resolution of HST spectra could also 
be used to obtain separate spectra of both the A and B 
sources to conclusively establish the origin of the three 
components. Likewise, spatially resolved spectra of the 
Paschen lines in the infrared (using adaptive optics from 
the ground, for instance) should be fruitful. In either 
case, we expect such observations to reveal that either A 
or B (but not both) would be associated with the broad 
emission-line profiles. 
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